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ABSTRACT 

HD 32297 is a young A-star (~30 Myr) 112 pc away with a bright edge-on debris disk that has been resolved in 
scattered light. We observed the HD 32297 debris disk in the far-infrared and sub-millimeter with the Herschel 
Space Observatory PACS and SPIRE instruments, populating the spectral energy distribution (SED) from 63 to 
500 jim. We aimed to determine the composition of dust grains in the HD 32297 disk through SED modeling, using 
geometrical constraints from the resolved imaging to break the degeneracies inherent in SED modeling. We found 
the best fitting SED model has two components: an outer ring centered around 1 10 AU, seen in the scattered light 
images, and an inner disk near the habitable zone of the star. The outer disk appears to be composed of grains >2 pm 
consisting of silicates, carbonaceous material, and water ice with an abundance ratio of 1:2:3 respectively and 90% 
porosity. These grains appear consistent with cometary grains, implying the underlying planetesimal population is 
dominated by comet-like bodies. We also discuss the 3.1a detection of [C ii] emission at 158 pm with the Herschel 
PACS instrument, making HD 32297 one of only a handful of debris disks with circumstellar gas detected. 
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1. INTRODUCTION 

Debris disks are circumstellar disks composed of dust pro- 
duced during collisions of planetesimals. In the youngest disks 
(10-100 Myr), the planetesimals may deliver volatiles such as 
water to still-forming terrestrial planets. Although the planetesi- 
mals themselves are undetectable, the dust they produce radiates 
thermally in the infrared (IR) and sub-millimeter (mm) and scat- 
ters starlight at optical and near-IR wavelengths. These grains 
provide clues to the composition of their parent bodies. 

The outer regions of several debris disks have been imaged in 
scattered light. Resolved images can constrain the morphology 
of the disk’s outer regions, but the composition of their grains 
cannot be uniquely determined from these images alone. Mid-IR 
spectra — useful for determining grain composition in younger 
protoplanetary disks — also fail to provide constraints on the dust 
composition in most debris disks; by this point, the remaining 
grains are too large to emit solid state features in the mid-IR 
(Chen et al. 2006). While there are a few notable exceptions, 
most debris disks require modeling of the full spectrum, from 
optical to mm wavelengths, to probe the grain composition. 

The Herschel Space Observatory (Pilbratt et al. 2010) was 
launched in 2009 May, presenting a new opportunity for 
sensitive far-IR and sub-mm observations. The PACS (Poglitsch 
et al. 2010) and SPIRE (Griffin et al. 2010) instruments 
have photometric and spectroscopic capabilities spanning a 
wavelength range of ~60-500 pm. These data are crucial for 
detailed modeling of a disk’s spectral energy distribution (SED) 
because they span the wavelength range where the thermal 
emission from debris disks typically peaks. 

Unfortunately, SED modeling of debris disks is hampered by 
degeneracies in the models between disk geometry and grain 
properties. Thankfully, resolved imagery can be used to break 
some of the degeneracies by providing geometrical constraints. 
We have used the Herschel Space Observatory's PACS and 


SPIRE instruments to obtain far-IR and sub-mm photometry 
and spectroscopy of the disk surrounding the ~30 Myr old 
(Kalas 2005) A-star, HD 32297. This edge-on disk has been 
imaged several times in scattered light, thereby constraining the 
disk geometry. The Herschel observations fill in a large gap in 
the SED, which allows us to model the grain composition in 
more detail. 

The HD 32297 disk (112 pc away; van Leeuwen 2007) was 
first resolved in the near-IR out to a distance of 3" 3 (400 AU) 
from the star with Hubble Space Telescope ( HST ) NICMOS 
(Schneider et al. 2005) and later resolved at several other 
near-IR wavelengths (Debes et al. 2009; Mawet et al. 2009). 
Recently, angular differential imaging (ADI) has been used 
to resolve the disk in the near-IR with ground-based facilities 
(Currie et al. 2012; Boccaletti et al. 2012). Additionally, the disk 
has been marginally resolved at mid-IR (Moerchen et al. 2007; 
Fitzgerald et al. 2007) and millimeter wavelengths (Maness et al. 
2008). 

The HD 32297 debris disk has a few unique features; one of 
the more luminous debris disks (L IR /L* ~ 1 0 3 ) , HD 32297 is 
also one of only a handful of debris disks where circumstellar gas 
has been detected. Redfield (2007) detected Na i in absorption 
toward HD 32297 that was not found toward any neighboring 
stars. The additional peculiarity of brightness asymmetries and 
warping seen in scattered light images were analyzed by Debes 
et al. (2009), who concluded that these features could be caused 
by the disk’s motion through the interstellar medium (ISM). 

The Herschel data presented here were acquired as part of 
the Herschel Open Time Key Programme entitled “Gas in Pro- 
toplanetary Systems” (GASPS; Dent et al. 2013). The PACS 
data were taken as part of the main program, and the SPIRE 
data were taken as part of an Open Time proposal to follow up 
GASPS debris disks at longer wavelengths (OT2_aroberge_3; 
PI: A. Roberge). Here we present the PACS and SPIRE obser- 
vations of HD 32297 as well as the results of the modeling of 
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the entire SED. In Section 2, we present the data and describe 
the data reduction. In Section 3, we describe the SED modeling 
and we discuss the results in Sections 4 and 5. 

2. OBSERVATIONS AND DATA REDUCTION 

Observations of HD 32297 at 70, 100, and 160 /xm were taken 
with the Herschel PACS instrument in ScanMap mode with the 
medium scan speed of 20" s _1 . The 70 /xm observations were 
taken at a scan angle of 63° and consisted of 8 scan legs with 
scan lengths of 3', and 2" separation between the legs, for a total 
observing time of220s.The 100 and 1 60 /im observations were 
taken simultaneously with slightly larger maps of 10 legs with 
the same leg length and separation as the 70 /im observation. 
Data at 100 and 160 /im were taken at two different scan angles, 
70 and 110°, with a total duration of 276 s per scan angle. 
The observations at the two scan angles were then combined to 
reduce noise due to streaking in the scan direction. The data were 
reduced with HIPE 8.2 (Ott 2010) using the standard reduction 
pipeline. The final maps were chosen to have a pixel scale 
corresponding to the native pixel scale of the PACS detectors, 
3"2 for the 70 and 100 /xm images and 6" 4 for the 160 /rm 
images. 

Simultaneous images at 250, 350, and 500 /xm of HD 32297 
were taken with the SPIRE instrument. The observations were 
made in the Small Scan Map mode with a scan speed of 30" s -1 , 
with two repetitions and a total observation time of 307 s. The 
data were reduced with HIPE 8.2, producing final maps with 
pixels scales of 6, 10, and 14" for the 250, 350, and 500 /xm 
images respectively. HIPE produces SPIRE images with units 
of Jy beam 1 , which we converted into Jy pixel 1 for analysis 
using beam areas of 423, 751, and 1587 arcsec 2 for the 250, 
350, and 500 /xm images respectively. 

The Herschel PACS spectroscopy of HD 32297 was taken in 
two modes, the lineSpec and rangeSpec modes. The lineSpec 
observations targeted the [Oi] 63.2 /xm line with a duration of 
3316 s. The rangeSpec observations targeted six lines, [Oi] 
at 145.5 /urn, [Cii] at 157.7 /xm, two 0 -H 2 O lines at 78.7 and 
179.5 /xm, and 2 CO lines at 72.8 and 90.2 /xm. The total 
rangeSpec observation time was 5141 s, divided into three 
observing segments, covering the six lines two at a time. A 
deep follow-up rangeSpec observation was performed targeting 
just the [Cii] 157.7 (im line with a duration of 4380 s. The 
spectroscopic data were reduced using HIPE 8.2. An aperture 
correction is applied in HIPE to account for point source flux 
loss. The spectra were produced with a pixel scale corresponding 
to the native resolution of the instrument. 

3. ANALYSIS 

3.1. Herschel PACS Photometry 

The HD 32297 disk is a spatially unresolved point source 
to Herschel. The full-width at half-maximum of the Herschel 
PACS beam at 70 /xm is 5"6, so the bulk of the thermal emission 
must be within ~300 AU. This is consistent with resolved 
images which suggest the disk peaks at ~110 AU (Debes 
et al. 2009; Currie et al. 2012; Boccaletti et al. 2012). Aperture 
photometry was performed with apertures of 12" for the 70 and 
100 /xm images and 22" for the 160 /xm image. 

Aperture corrections provided by the Herschel PACS ICC 5 
were applied, but color corrections were not applied. The 

5 PICC-ME-TN-037: http://herschel.esac.esa.int/twiki/pub/Public/ 
PacsCalibrationWeb/pacs_bolo_fluxcal_report_vl.pdf 


Table 1 

Herschel PACS and SPIRE Photometry Results 


ObsID 

Wavelength 
O m) 

Flux ± Error 

(Jy) 

1342193125 

70 

1.038 ±0.029 

1342217452-3 

100 

0.770 ± 0.022 

1342217452-3 

160 

0.403 ± 0.020 

1342240033 

250 

0.153 ±0.012 

1342240033 

350 

0.071 ±0.008 

1342240033 

500 

0.045 ± 0.007 


uncertainty in the flux was calculated from the standard devia- 
tion of the sky background in an annulus around the aperture. 
The annulus was placed between 20" and 30" from the central 
star for the 70 and 100 /xm images and between 30" and 40" 
for the 160 /xm image. An absolute calibration error of 2.64%, 
2.75%, and 4.15% (see footnote 5) for the 70, 100, and 160 /xm 
images respectively was added in quadrature to the uncertainty 
measured from the sky background. Results from the Herschel 
PACS photometry observations are given in Table 1 . 

3.2. Herschel SPIRE Photometry 

Aperture photometry was performed with aperture radii of 22, 
30, and 42" for the 250, 350, and 500 /xm images respectively. 
The sky background was estimated from an annulus with 
radius of 60-90" from the central star and subtracted from the 
measured flux. Aperture corrections were applied according 
to the SPIRE data reduction guide. 5 A color correction of 
~5% was applied assuming a Rayleigh-Jeans law slope of 
F v ex v 2 . This correction could be off by 5% if the slope 
of the SED in the sub-mm is steeper than Rayleigh-Jeans. The 
uncertainties in the flux measurements come from the standard 
deviation of the sky background added in quadrature with a 
7% absolute calibration error. Results from the Herschel SPIRE 
photometry are listed in Table 1 . 

3.3. Herschel PACS Spectroscopy 

The PACS instrument is an Integral Field Unit (IFU) spec- 
trometer that has a 5 x 5 array of spaxels with each spaxel cover- 
ing a 9'.' 4 x 9'.' 4 region. We verified that the star was well-centered 
on the central spaxel during the observations by comparing the 
observations to a model of the transmission of a theoretical 
point-spread function (PSF) through the PACS IFU. We shifted 
the model PSF and calculated the fractional flux in the differ- 
ent spaxels as a function of offset. Comparison to observations 
of HD 32297 indicate the star is not significantly offset. This 
method is the same as the one used by C. D. Howard et al. 
(2013, in preparation) for PACS observations of Taurus that had 
pointing errors. 

We use data from the central spaxel only for analysis. No lines 
were detected in the lineSpec or the first rangeSpec observation. 
We calculated continuum flux values by fitting straight lines 
to the data, and the results are listed in Table 2. We excluded 
five pixels on each edge from the continuum fit because PACS 
spectra have enhanced noise at the edges. Emission in the 
[Cii] 157.7 /xm line was seen in the second deep rangeSpec 
observation. The continuum flux was found by fitting a line to 
the spectrum while masking out 0.5 /xm around the line center. 
The rms noise was also calculated in the region surrounding 
the line. The number of pixels used to calculate the noise was 

6 SDRG5.7: http://herschel.esac.esa.int/hcss-doc-9.0/ 
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Figure 1. The [C n] 157.7 /im emission line from the HD 32297 debris disk. The 
dashed line represents the expected position of the [C n] line. The continuum 
was fit with a straight line (black solid line), and the emission line was fit with 
a Gaussian (red solid line). The integrated line flux was measured from the 
three pixels surrounding the line (indicated by the gray bar). The integrated line 
flux is 3.7cr above the continuum. 

(A color version of this figure is available in the online journal.) 


Table 2 

Herschel PACS Spectroscopy Results 


Line 

Name 

Wavelength 

(/xm) 

Continuum 
Flux (Jy) 

Integrated Line Flux 3 
x 10~ 18 (Wm~ 2 ) 

Oi 63 

63.184 

1.28 ±0.20 

<7.29 

CO 72 

72.843 

0.83 ± 0.23 

<7.77 

H 2 0 79 

78.741 

1.03 ±0.23 

<9.12 

CO 90 

90.163 

1.05 ±0.29 

<8.13 

Oi 145 

145.525 

0.55 ± 0.09 

<3.99 

Cii 158(l) b 

157.741 

0.53 ±0.13 

<4.32 

H 2 0 180 

179.741 

0.44 ±0.15 

<3.57 

Cii 158 (2) b 

157.741 

0.50 ± 0.06 

2.68 ±0.72 


Notes. 

3 Non-detections are reported as 3<r upper limits. 

b Cn 158 line was observed twice: (1) first rangeSpec observation (2) second 
deeper rangeSpec observation. 



T ex (K) 

Figure 2. Column density of Cn as a function of excitation temperature 
(solid line). The dashed line shows the lower limit on the column density 
( N = 2.5 x 10 11 cm~ 2 ) for excitation temperatures <300 K. 

where the indices 1 and 0 indicate the upper and lower fine- 
structure levels respectively, X = 157.7 /tm, F io = 2.68 x 
10~ 18 W m~ 2 is the integrated line flux, A,o = 2.4 x 10~ 6 s -1 
is the spontaneous transition probability, £2 = 0.357 arcsec 2 is 
the angular source size estimated from resolved scattered light 
images, and x\ is the fractional population of the upper level. 
Assuming local thermodynamic equilibrium, v, is dependent on 
the excitation temperature. Roberge & Kamp (2011) give this 
as 

(27, + 1 )g x e- E '/ kT - 

x\ = , (2) 

Q(T ex ) 

where 7, = 3/2 is the angular momentum quantum number 
of the upper level, g.v are the nuclear statistical weights for 
the two levels (go = 2 and g, = 4), E\/k = 91.21 K is the 
energy difference between the two levels, T ex is the excitation 
temperature, and Q(T ex ) is the partition function, which here 
can be approximated in a two-level system by 

Q(T ex ) = go + g l e- E1/kn \ (3) 


chosen such that the signal-to-noise in the emission line was 
maximized, i.e., the rms noise was minimized. 

The line flux was integrated over the three pixels surrounding 
the line (marked by the gray bar in Figure 1). The uncertainty 
in the line flux was found by propagating the rms noise. Upper 
limits on the other lines were found using the same method of 
error propagation at the expected line center. These values are 
reported in Table 2. The integrated line flux of the [C ii] line is 
detected at the 3.7 a level above the continuum in the central 
spaxel from the deep rangeSpec observation. All 25 spaxels 
were searched for a significant [C ii] line signal, and no other 
emission was found with more than a 2 a significance. Since the 
line is present only in the central spaxel, we believe the gas is 
associated with the star rather than from the surrounding ISM. 

3.4. Column Density ofCn in HD 32297 

With the flux of the [Cii] line, we can determine the 
column density of Cii for comparison to the Nai column 
density found by Redfield (2007). The [C ii] 158 pm line arises 
from the transition between the two fine structure lines of the 
electronic ground state. Since the next electronic energy level 
is three orders of magnitude higher, we can assume a two-level 
population with a high accuracy. Following Roberge & Kamp 
(2011), the column density is 


he A,ox,£2 


With only one gas line, we cannot measure the excitation 
temperature. However, we can derive a lower limit on the column 
density. Figure 2 shows the dependence of the column density 
on excitation temperature. From 1 to 300 K, we find a lower 
limit on the Cii column density of Nc u > 2.5 x 10 11 cm -2 . 
This value is similar to the column density of Na i calculated by 
Redfield (2007) (JV NaI = 2.5 x 10 11 cm" 2 ). 

4. SED MODELING 
4.1. SED Data 

In addition to the Herschel data, we collected archive data 
for use in our SED modeling (see Table 3). To constrain the 
stellar photosphere, we used Hipparcos B and V (Perryman & 
ESA 1997), 2MASS 7, //, and K s (Cutri et al. 2003), and WISE 
bands 1 and 2 (Wright et al. 2010). For the IR excess, we used 
IRAS 25 and 60 pm (Moor et al. 2006), WISE bands 3 and 4 
(Wright et al. 2010), Spitzer MIPS 24 pm (Maness et al. 2008), 
and Spitzer IRS data taken from the enhanced data product in 
the Spitzer Heritage Archive. 7 We also included millimeter data 
from Meeus et al. (2012): the 1.2 mm flux from the MAMB02 
bolometer array on the IRAM 30 m telescope and the 1.3 mm 
flux from the Sub-Millimeter Array (SMA). The uncertainties 
for these last two measurements (reported in Table 3) include 
15% calibration uncertainties added in quadrature. 


irsa. ipac . caltech . edu/ applications/Spitzer/SH A/ 
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Table 3 

Additional Data Used in SED Modeling 


Wavelength 

(^m) 

Flux and Uncertainty 

(Jy) 

Instrument 

Reference 

0.438 

1 .952 ± 0.026 

Hipparcos B 

i 

0.547 

2.094 ± 0.023 

Hipparcos V 

i 

1.235 

1 .342 ± 0.030 

2MASS J 

2 

1.65 

0.913 ± 0.044 

2MASS H 

2 

2.16 

0.611 ±0.010 

2MASS K s 

2 

3.4 

0.278 ±0.010 

WISE 1 

3 

4.6 

0.150 ±0.005 

WISE 2 

3 

8.00 

0.064 ± 0.002 

Spitzer IRS 

4 

8.65 

0.059 ± 0.003 

Spitzer IRS 

4 

9.35 

0.053 ± 0.003 

Spitzer IRS 

4 

10.11 

0.049 ± 0.002 

Spitzer IRS 

4 

10.93 

0.049 ± 0.002 

Spitzer IRS 

4 

11.2 

0.050 ± 0.002 

Gemini-N /Michelle 

5 

11.56 

0.053 ± 0.005 

WISE 3 

3 

11.7 

0.053 ± 0.005 

Gemini-S /T-ReCS 

6 

11.81 

0.050 ± 0.002 

Spitzer IRS 

4 

12.77 

0.052 ± 0.001 

Spitzer IRS 

4 

13.80 

0.054 ± 0.003 

Spitzer IRS 

4 

14.92 

0.063 ± 0.002 

Spitzer IRS 

4 

16.13 

0.072 ± 0.003 

Spitzer IRS 

4 

17.44 

0.087 ± 0.004 

Spitzer IRS 

4 

18.3 

0.090 ±0.014 

Gemini-S /T-ReCS 

6 

18.85 

0.110 ±0.004 

Spitzer IRS 

4 

20.38 

0.142 ±0.009 

Spitzer IRS 

4 

22 

0.193 ±0.020 

WISE A 

3 

22.03 

0.189 ±0.009 

Spitzer IRS 

4 

23.68 

0.210 ±0.010 

Spitzer MIPS 

7 

23.81 

0.239 ±0.011 

Spitzer IRS 

4 

25 

0.256 ± 0.041 

IRAS 25 

8 

25.74 

0.296 ±0.010 

Spitzer IRS 

4 

27.83 

0.375 ±0.011 

Spitzer IRS 

4 

30.08 

0.444 ±0.011 

Spitzer IRS 

4 

32.51 

0.563 ± 0.014 

Spitzer IRS 

4 

35.15 

0.668 ±0.018 

Spitzer IRS 

4 

60 

1.140 ±0.070 

IRAS 60 

8 

1200 

0.00314 ±0.00095 

IRAM 30m/MAMBO2 

9 

1300 

0.00310 ±0.00087 

SMA 

9 


Notes. Color-corrected flux with la error bars used in the SED modeling. 
References. ( 1) H0g et al. 2000; (2) Cutri et al. 2003; (3) Wright et al. 2010; (4) 
Spitzer Heritage Archive; (5) Fitzgerald et al. 2007; (6) Moerchen et al. 2007; 
(7) Maness et al. 2008; (8) Moor et al. 2006; (9) Meeus et al. 2012. 

The HD 32297 disk is also marginally resolved at mid-IR and 
millimeter wavelengths. We used the total flux from Gemini 
North Michelle imaging at 1 1 .2 /xm (Fitzgerald et al. 2007), and 
Gemini South T-ReCS imaging at 1 1.7 and 18.3 /xm (Moerchen 
et al. 2007). We used the SMA flux at 1.3 mm rather than the 
total flux from the CARMA 1.3 mm resolved image (Maness 
et al. 2008), because the unresolved data from the SMA has 
a smaller uncertainty. The photometric data used are listed in 
Table 3 and plotted in Figure 3. 

We combined these data with our Herschel PACS and 
SPIRE photometry and the continuum values from the PACS 
spectroscopy. Our data points fill in a large gap in the SED from 
60 to 500 /xm (see Figure 3) and allow us to assess where the 
peak of the thermal emission is located. 

4.2. Stellar Properties 

We started our analysis by fitting for the stellar parameters 
of HD 32297 . We used the Hipparcos B and V, 2MASS J, H, 
and K s , and WISE bands 1 and 2 to constrain the stellar models. 



Figure 3. Spectral energy distribution of HD 32297 with simple blackbody fits. 
The gray line is the best-fitting stellar photosphere model with T = 7750 K 
(see Section 4.2). The data points that were taken from the literature and listed 
in Table 3 are plotted here as well as the new PACS and SPIRE data given in 
Tables 1 and 2. A two-temperature modified blackbody model is overplotted 
(orange solid line), with the individual components also shown — a wanner 
blackbody with T = 240 K (dotted line) and a colder blackbody with T = 83 K 
(dashed line) as described in Section 4.4. 

We fit the stellar data with ATLAS9 stellar photosphere models 
(Castelli & Kurucz 2004). The spectral type of HD 32297 is 
usually quoted as either an A0 (Torres et al. 2006) or an A5 
(Fitzgerald et al. 2007, and references therein). Our best fitting 
model with an A0V spectral type ( T — 9750 K) is too hot and 
does not match the data well. 

We therefore tried both adding interstellar extinction with 
a Fitzpatrick (1999) extinction law and varying the stellar 
photosphere temperature. Unfortunately, these two parameters 
are degenerate. Fixing the temperature to T = 9750 K, the 
best fitting model has an extinction of Ay = 0.63 ± 0.02 mag, 
a reduced chi-squared value of xf — 0.59, and a bolometric 
luminosity of L = 11.9 L 0 . The best fitting model with ex- 
tinction and temperature as free parameters has T = 8000 K, 
A v = 0.161 ± 0.026 mag, and xl = 0.55. 

To break the degeneracy, we added UV continuum values at 
0.26 and 0.31 /xm from unpublished STIS spectra (S. Redfield 
et al. 2012, in preparation). The best model is one with no ex- 
tinction, a temperature of T = 7750 K, and a luminosity of L = 
5.6 L q (Figure 3). This is more consistent with an A7 spectral 
type than an A0. This hint in the UV spectra that the spectral type 
is later than was thought will be thoroughly analyzed in a future 
paper (S. Redfield et al. 2012, in preparation). We note that this 
lower temperature is the same used by Debes et al. (2009) and 
close to the temperature used by Fitzgerald et al. (2007), both 
of whom note that the star appears to be under-luminous for a 
main sequence star of this temperature. Some reasons proposed 
are errors in the Hipparcos distance (Debes et al. 2009) or a 
sub-solar metallicity (Meeus et al. 2012). 

The assumed stellar luminosity does not have a significant 
effect on the SED in the far-IR where the stellar contribution 
is negligible. However, the assumed temperature of the star is 
important for calculating the equilibrium temperature of the 
dust. An error in the stellar temperature will translate into an 
error in the dust temperature, which will affect the calculated 
grain properties, such as minimum grain size. 

4.3. Surface Density Profile 

The HD 32297 disk has been well resolved in scattered light 
with 775TNICMOS at 1.1 /xm (Schneider et al. 2005), 1.6 /xm, 
2.05 /xm (Debes et al. 2009), and K s band (7. = 2.16 /xm) from 
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the ground with Keck/NIRC2 (Currie et al. 2012), as well as 
H (X = 1.65 /xm) and K s bands with Very Large Telescope/ 
NACO (Boccaletti et al. 2012). These observations block out 
the central star with a coronagraph, and consequently obscure 
the inner portions of the disk as well. 

The resolved images place strong geometrical constraints 
on the disk outside of ~65 AU. Currie et al. (2012) and 
Boccaletti et al. (2012) have both recently published models 
of the disk based on their near-IR ground based imaging using 
ADI. Boccaletti et al. (2012) warn that images processed using 
ADI and/or the Locally Optimized Combination of Images 
(LOCI) are subject to self-subtraction and other artifacts (also 
see Milli et al. 2012). Therefore, a direct inversion of the surface 
brightness profile in the manner of Augereau & Beust (2006) is 
not practical here. For this reason, we rely on the disk modeling 
that takes into account the ADI and LOCI processing. 

The models of Currie et al. (2012) and Boccaletti et al. 
(2012) agree quite well. Both see a break in the surface 
brightness profiles around 1" (~110 AU), which was also seen 
by Debes et al. (2009) in the NICMOS images. Both measure an 
inclination of 88°, i.e., 2° from edge-on. This is notably different 
from the inclination measured by Schneider et al. (2005) of 79?5. 

Where the models disagree is in the anisotropic scatter- 
ing factor, g. Currie et al. (2012) use a two-component 
Henyey-Greenstein phase function with a highly forward scat- 
tering component (gi = 0.96) and a backscattering compo- 
nent (g 2 = —0.1), while Boccaletti et al. (2012) use only a 
one-component phase function with a best fit value of g — 0.5. 
They discuss how higher values of g would make the disk too 
bright in the inner regions. Currie et al. (2012) also note this, but 
they dismiss it due to the large uncertainties in surface brightness 
close to the star. 

We chose to use the models of Boccaletti et al. (2012) because 
they used the GRaTer code (Augereau et al. 1999; Lebreton 
et al. 2012) to model the disk images, which we also used 
to model the SED (see Section 4.5). The best-fit model of 
Boccaletti et al. (2012) to the K s band image has a mid-plane 
density of the form 

„m = ».V2 'TV' 1 . (« 

V ' \V 110 AU/ V 110 AU/ ) 

We assumed the disk is geometrically thin because the SED 
modeling cannot distinguish between a vertical offset from the 
midplane and a radial change in distance. We chose to use the 
Boccaletti et al. (2012) model of the K s band image rather than 
the H band image because the K s band data are of better quality. 

4.4. Dust Disk Modeling Strategy 

We started our modeling by fitting the IR excess with a single- 
temperature modified blackbody, i.e., a blackbody model with 
an opacity at longer wavelengths of the form iV . The exact form 
of the modified blackbody we used is F v oc r v B v for A. 7 ,q, and 
F v oc B v for /, < 7,(|, where B v is a simple blackbody and r is the 
optical depth, which takes the form r oc ( v/vq )P oc (X/Xq)~^ . 
Here we have fixed [5 = 1 and Xo = 1 00 /iin, based on typical 
values for debris disks (Dent et al. 2000; Williams & Andrews 
2006). The best fit was determined through x 2 minimization 
using MPFIT (Markwardt 2009). The best fitting model has a 
temperature of T = 80 K. However, this model significantly 
underestimates the flux at mid-IR wavelengths. 

To account for the missing IR flux in our model, we added 
a second modified blackbody. The best fitting model has two 


Table 4 

Parameters Explored in GRaTer Models 


Parameter 

Range 

No. of Values 

Distribution 

K 

-5 to -2.5 

20 

Linear 

a min 

0.01 to 100 u m 

77 

Logarithmic 

Carbon volume 

0% to 100% 

21 

Linear 

Ice volume 

0% to 90% 

10 

Linear 

Porosity 

0% to 95% 

20 

Linear 


distinct components with temperatures of T = 83 K and 
T — 240 K (see Figure 3). The addition of the second blackbody 
much improves the fit over the single-temperature model, 
suggesting there is a second inner disk. The inner component is 
too hot, and therefore too close to the star, to be part of the outer 
ring imaged in scattered light. 

Spatial information is important for breaking degeneracies in 
SED modeling. Scattered light images of HD 32297 restrict 
the models of the outer disk geometry. Unfortunately, no 
constraining images exist for the inner disk, which is too close 
to the star to be imaged. Therefore, we divided the dust disk 
modeling into two steps: first fitting the outer disk with a more 
complex model, then fitting the poorly constrained inner disk 
with a simpler one. 

4.5. Outer Disk Modeling with GRaTer 

For the outer disk, we limited the data to those with wave- 
lengths larger than 25 /urn. This cutoff was chosen as the point 
where the inner disk contributes less than 50% to the total 
two-component blackbody model. With the Herschel observa- 
tions, we have enough data points with wavelengths greater than 
25 /rm to constrain the outer disk properties. 

We use the GRaTer code (Augereau et al. 1999; Lebreton 
et al. 2012) to fit the SED of HD 32297. GRaTer is a dust 
disk modeling code specifically designed for optically thin 
debris disks, which can compute large grids of models with 
different grain sizes and composition and use a density profile 
to describe the disk geometry. The code computes both the 
scattered light emission and the thermal emission from the grains 
in equilibrium with the radiation field of the central star. We 
constrained the spatial distribution of the dust grains by using 
the models of the resolved scattered light images from Boccaletti 
et al. (2012). We confined the grains in our SED model to be 
located in a ring defined by Equation (4). The overall abundance 
was left as a free parameter that scales the radial profile to match 
the SED. 

With the disk geometry fixed, we focused our modeling on the 
grain sizes and composition, which we assumed to be the same 
throughout the disk. GRaTer calculates Mie scattering and 
absorption coefficients for a large range of grain compositions. 
Specifically, we explored combinations of materials consisting 
of astrosilicate (Draine 2003), amorphous carbon (Zubko et al. 
1 996), amorphous water ice (Li & Greenberg 1 998) and porosity. 
The volume ratios of the materials explored are listed in Table 4. 
If any material reaches its sublimation temperature, it is replaced 
by vacuum. For most of the outer disk, the temperatures are too 
cold for this to happen. 

For the grain sizes, we explored a range of minimum grain 
sizes from a m ; n = 0.01 /xm to a m = 100 fim and grain size 
distributions with a power-law of the form n(a)da oc cT da, with 
k ranging from —5 to —2.5. The maximum grain size is kept 
fixed at a max = 7.8 mm, which is large enough compared to the 
longest wavelength of the data (X = 1.3 mm). However, we only 
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Figure 4. Top: spectral energy distribution for the best-fitting total model given in Table 5 (Model 6). The black solid line is the fit to the stellar photosphere, the 
dotted and dashed lines show the inner and outer disk models respectively, and the red solid line represents the combined model of photosphere + inner disk + 
outer disk. The data plotted here are listed in Tables 1—3. Bottom: the relative residuals of the best-fitting model shown above. The relative residuals are of the form 
(data - model#6)/data. Also plotted for comparison are the relative residuals of the other five models given in Table 5, e.g., (model#l - model#6)/model#l. The red 
line marks the zero point, or the residual for model 6. 


Table 5 

Results of SED Modeling 


Outer Disk Composition 3 

Volume 

Ratios 

Density 
(g cm" 3 ) 

Outer Disk 

Inner Disk 


Total Disk 

^min 

K 

r min 

*2 min 

xl 

V 

Dust Mass b 

Model 1 (AS) 


3.50 

1.4 ± 0.1 /zm 

-3.8 ±0.2 

3.3 ± 0.3 AU 

200 ± 120 /im 

4.58 

36 

0.32 ± 0.05 M e 

Model 2 (AS+C) 

1:4 

2.46 

2.1 ± 0.4 /zm 

-4.3 ±0.3 

3.5 ± 0.5 AU 

31.6 ± 26 /xm 

4.27 

35 

0.11 ±0.02 

Model 3 (AS+P) 

11:9 

1.93 

0.3 ±0.1 /xm 

-3.5 ±0.2 

3.6 ± 0.4 AU 

31.6 ± 27 /xm 

4.42 

35 

0.26 ± 0.05 M e 

Model 4 (AS+I) 

4:1 

3.04 

1.0 ± 0.1 /xm 

-3.7 ±0.2 

3.2 ± 0.3 AU 

39.8 ±22 /zm 

4.45 

35 

0.34 ± 0.12 M e 

Model 5 (AS+C+P) 

1:2:12 

0.53 

3.4 ± 0.1 /zm 

-3.8 ±0.2 

2.8 ± 0.3 AU 

25.1 ± 24 /um 

3.60 

35 

0.08 ± 0.01 M e 

Model 6 (AS+C+I+P) 

1:2:3:54 

0.19 

2.1 ± 0.3 /zm 

-3.3 ±0.2 

1.1 ±0.2 AU 

2.2 ±0.9 Atm 

1.59 

34 

0.10 ± 0.01 M e 


Notes. 

a Model 1 : Astrosilicate (AS), Model 2: Astrosilicate + Carbon (C), Model 3: Astrosilicate + Vacuum (P), Model 4: Astrosilicate + Water Ice (I), Model 5: Astrosilicate + 
Carbon + Vacuum, Model 6: Astrosilicate + Carbon + Water Ice + Vacuum. 

b The dust mass is calculated for dust grains smaller than 1 mm only. The total dust mass is dominated by the outer disk. 


considered grains with sizes smaller than 1 mm in calculating 
the dust mass (given in Table 5), to be consistent with other 
results from GASPS and other Herschel Key Programmes. 

We modeled the disk using six different compositions with 
varying levels of complexity. The simplest model used only 
astrosilicate grains. The next three models used silicates mixed 
with only one other grain type: carbonaceous material, water 
ice, or increased porosity. For the three and four material 
combinations, we kept the silicate to carbon volume ratio fixed 
at 1:2. This is the ratio expected from cosmic abundances and 
is similar to the ratio observed in comet Flalley dust (Greenberg 
1998). 

The grains are assumed to be porous aggregates of silicate, 
carbon and water ice. The scattering and absorption coefficients 
for the aggregates are calculated using the Bruggeman mixing 
rule. Silicate and carbon are mixed first, then the Si+C mixture is 
mixed with water ice, and finally it is all mixed with vacuum to 
simulate porous grains. For more details see Lebreton et al. 
(2012) and Augereau et al. (1999). Grain densities for the 
mixtures appear in Table 5. 


4.6. Inner Disk Modeling 

The mid-IR flux can only be explained with a warmer 
component than is seen in the scattered light images. We first 
tried to model the disk with the geometry constrained only by 
the coronagraphic scattered light images of the disk. We found 
that the thermal emission from the grains seen in the scattered 
light images was not enough to reproduce the mid-IR flux in 
the SED. The data in this region come from several sources, 
and all agree with each other ( Spitzer IRS, WISE (Wright et al. 
2010), and Gemini N (Fitzgerald et al. 2007) and S (Moerchen 
et al. 2007)). Therefore, we added a warmer component to fit 
the mid-IR data. 

We chose to model the warm component as an inner disk 
inside the radius masked by coronagraphs in the images. The 
inner disk is less constrained than the outer disk since the inner 
disk lacks geometrical information and the Lir/L* of the inner 
disk is an order of magnitude lower than the outer disk (6.9 x 
1CT 4 versus 5.6 x 10~ 3 ; see Figure 4). This leads to a degeneracy 
between two of the dust disk modeling parameters: the minimum 
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radius and the minimum grain size. The degeneracy exists 
because decreasing both the parameters has the same effect of 
increasing the grain temperature. Changing the minimum grain 
size also affects the amount of flux emitted from a given grain, 
but this can be mimicked by varying the total number of grains. 
With a weakly-emitting disk, the separate effects can be difficult 
to disentangle. 

For each outer disk model with a different grain composition, 
we subtracted it from the SED and fit our inner disk model to 
the residuals. Since we lack constraining spatial information for 
the inner disk, we needed to use a simpler model. We used the 
model described in Donaldson et al. (2012), which calculates 
only the thermal emission from astrosilicate grains in radiative 
equilibrium with the central star. We fixed the outer radius to 
5 AU since all values above this had no significant effect on the 
SED. We assumed the disk has a surface density profile of the 
form Xtr) oc r _1 5 , consistent with collision-dominated disks 
(Krivov et al. 2006; Strubbe & Chiang 2006). We also fixed 
the grain size distribution throughout the disk to a Dohnanyi 
(1969) power-law (ri(a)da oc a~ 35 da) with a maximum grain 
size of 1 mm. The only free parameters were the inner radius, 
the minimum grain size, and the dust mass. 

We used the inner disk model to fit the residuals after 
subtraction of the outer disk model, then combined this model 
with those of the outer disk and the star. We calculated reduced 
y 2 values (y 2 ) using all 43 data points listed in Tables 1-3 
greater than 8 p m and the number of degrees of freedom (v) 
listed in Table 5. We calculated the errors on the parameters 
from the ler confidence intervals in the y 2 distribution. 

A similar modeling approach was recently used by Ertel et al. 
(201 1) for HD 107146. They fit the disk SED using the scattered 
light images to constrain the geometry, and they also found 
an overabundance of Spitzer IRS mid-IR flux in their best fit 
models. They tried modeling the excess in two ways: first by 
adding a small grain population within the imaged disk, and 
second by adding an inner disk. The small grain model was 
unable to reproduce the mid-IR flux in HD 107146; the inner 
disk model was needed to match the Spitzer IRS spectrum. 

4.7. Results 

We found the best fitting outer disk model was the four- 
material composition of silicates, carbon, and water ice in a 1:2:3 
ratio with a high porosity of 90% (final row of Table 5, Figure 4). 
The outer disk grains have a minimum size of c/ mln = 2. 1 pm, 
with a grain size distribution power-law index of k = —3.3. 
The best fitting total SED model also includes an inner disk 
from 1.1 AU with an unconstrained outer edge. The inner disk 
was fit with astrosilicate grains with a Dohnanyi (1969) size 
distribution ( k = —3.5). The minimum grain size in the inner 
disk for the best fitting model is a m j n = 2.2 [i m, similar to the 
outer disk grains. Other grain models tested appear in Table 5. 
We determined the uncertainties given in Table 5 from the ler 
confidence intervals in the y 2 distribution after fixing the other 
free parameters to the value that gives the smallest y 2 value. 

The four-material composition model has the best reduced 
y 2 value of the tested compositions (y 2 = 1.59, v = 34). 
Of the two-material composition models, astrosilicate + carbon 
and astrosilicate + porosity are the best fitting (y 2 = 4.27 and 
y 2 = 4.42, v = 35), but the model is much improved when all 
three materials are added together (y 2 = 3.60, v = 35). The 
addition of ice makes a dramatic improvement in the fit from 
the three-material composition to the four-material composition. 
The minimum grain size in the inner disk models decreases with 


each improvement on the outer disk model, moving closer to the 
expected blowout size of 1 pm. This is likely because the two- 
material models fail to fit all the mid-IR flux coming from the 
outer ring, and the inner disk models have to compensate. 

The blowout size given above is calculated for the inner 
disk only, where the assumed grain composition is astrosilicate. 
Given a density of p = 3.5 g cm -3 , the blowout size for spherical 
grains on circular orbits (in microns) is «biow =1.15 L*/(M*p), 
with L* and M* in solar units and p in g cm 3 . For astrosilicate 
grains around HD 32297 (M* = 1.84 M 0 ,L* = 5.3 L 0 ), the 
blowout size is «biow = 1 /um. 

For the outer disk, the grains are highly porous and likely have 
a fractal structure. The calculation of the blowout size depends 
on the surface area of the grains; this is non- trivial for fractal 
grains. The above equation is for spherical grains only; we do 
not calculate the blowout size for the outer disk because it is not 
very realistic in this case. For more about this problem, see the 
discussion in Lebreton et al. (2012). 

5. DISCUSSION 

The composition of dust grains in young debris disks is a 
key piece in understanding the last stages of terrestrial planet 
formation. A handful of debris disks, including ft Pictoris 
and HD 172555 (Telesco & Knacke 1991; Lisse'et al. 2009), 
have solid state features in the mid-IR that indicate dust grain 
composition. But most debris disks, including HD 32297, lack 
these features and, therefore, modeling of the full SED is needed 
to constrain the grain composition. 

Unfortunately, the presence of an unresolved warm compo- 
nent to the HD 32297 system complicates the modeling. Without 
resolved imaging of the inner regions, it is impossible to know 
the distribution of the warm component and how much mid-IR 
flux is coming from the outer disk versus the inner disk. The 
models of the inner disk depend strongly on the model chosen 
for the outer disk. Additionally, since there are no geometri- 
cal constraints on the inner disk, the results depend upon the 
distribution we assumed. 

Another concern is the age of the system. Kalas (2005) states 
the age as 30 Myr based on an uncertain association with either 
the Gould belt or Taurus Aurigae. The dust mass of HD 32297 
(see Table 5) is high for a debris disk. The mass is comparable to 
the 8 Myr old disk HR4796A, which suggests HD 32297 may be 
younger than 30 Myr. HD 32297 is also one of the oldest debris 
disks with gas detected, another indication it may be a younger 
system. But the system is likely not much younger than the given 
age of 30 Myr. There are several indicators that HD 32297 is a 
main sequence star with an optically thin debris disk, including 
(1) the lack of optical extinction, (2) the low fractional dust 
luminosity (L ir /L* ~ 10~ 3 ) and (3) a lack of solid state features 
in the mid-IR indicating no small grains are present. 

5.1. Cometary Dust? 

The best fitting model to the outer disk includes grains 
that are highly porous and icy. High grain porosity is seen in 
interplanetary dust particles collected with Stardust (Brownlee 
et al. 2006) and in the ejecta from comet Temple 1 created by 
Deep Impact (A’Hearn 2008). Greenberg & Hage (1990) showed 
that comet Halley’s spectrum could only be fit by highly porous 
grains with a porosity between 93% and 97.5%. Li & Greenberg 
(1998) modeled the ft Pictoris disk with similar composition 
dust. They assumed the ft Pic dust was cometary in origin and 
rejected models of compact grains with porosity lower than 90%. 
Polarized light observations of AU Mic also indicate that disk is 
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dominated by highly porous grains with porosity of 91%-94% 
(Graham et al. 2007; Fitzgerald et al. 2007). The highly porous, 
icy dust around HD 32297 is similar to ft Pic and solar system 
comets. 

Dust in the outer ring of HD 32297 therefore appears 
consistent with cometary dust particles. The ring is centered 
around ~110 AU, far from the star where ices should be 
prevalent. If planetesimal collisions produce the dust, this 
indicates that comet-like bodies dominate the planetesimal 
population in the outer disk. 

A large population of comets in the outer disk could deliver 
water to terrestrial planets. A significant fraction of Earth’s 
water likely came from Kuiper Belt comets (Morbidelli et al. 
2000). At an age of 30 Myr, HD 32297 could still be forming 
terrestrial planets (Kenyon & Bromley 2006). Comets scattered 
into the inner regions of the disk could deliver water to forming 
terrestrial planets in the habitable zone. 

5.2. Grain Porosity and ISM Interaction 

The scattered light images of HD 32297 , as well as those of 
the edge-on disks HD 15115 and HD 61005, have asymmetric 
features that may be due to interaction with ISM as the systems 
move with respect to their surroundings. The short wavelength 
of images of HD 32297 (Schneider et al. 2005; Debes et al. 2009) 
and HD 61005 (Hines et al. 2007; Maness et al. 2009; Buenzli 
et al. 2010) show a “swept-out” feature, while HD 15115 has 
a strong east-west asymmetry (Kalas et al. 2007; Debes et al. 
2008; Rodigas et al. 2012). The ISM interaction model of Debes 
et al. (2009) reproduces the features of all three disks. 

The ISM affects the disk grains through gas drag and/or 
grain-grain collision as the disk moves through a dense clump 
in the ISM. Unbound or weakly bound grains are swept back as 
they interact with the ISM. The most affected grains are typically 
thought to be the small grains that are nearly unbound due to 
the effect of radiation pressure. But larger grains with a higher 
porosity are also strongly affected by radiation pressure, and 
therefore may also be susceptible to being blown back through 
ISM ram pressure. The high porosity of the outer disk grains 
may be one factor that helps explain why this disk has such a 
strong ISM interaction feature. Stellar motion and environment 
must also play a role, since other disks modeled with a high 
grain porosity, such as /3 Pic, do not have the same feature. 

5.3. Gas in HD 32297 

HD 32297 is one of only a handful of debris disks that have 
gas detections. Redfield (2007) found Nai in absorption, aided 
by the disk’s nearly edge-on orientation. The detection of [C ii] 
emission from HD 32297 is the fourth detection of atomic gas 
from a debris disk with Herschel, though it is weaker than the 
lines seen from ft Pictoris, HD 172555, and 49 Ceti (Brandeker 
et al. 2012; Riviere-Marichalar et al. 2012; Roberge et al. 2013). 
It is also unusual that [C ii] was detected while [O i] was not. 
The only other debris disk with gas where this is true is 49 Ceti 
(Roberge et al. 2013). 

Given relatively advanced age of HD 32297 (~30 Myr; Kalas 
2005) and the typical protoplanetary disk lifetime (<10 Myr; 
Williams & Cieza 2011), the HD 32297 gas is unlikely to be 
primordial. The lack of sub-mm CO emission suggests a gas- 
to-dust ratio lower than is seen in younger protoplanetary disks 
(e.g., Zuckerman et al. 1995). Furthermore, the disk dust has a 
relatively low abundance, and shares other characteristics with 
debris dust (no detectable line-of-sight extinction, a lack of 


mid-IR solid state features from small grains). With little dust 
and little gas, the disk should be optically thin to stellar and 
interstellar dissociating U V photons and molecular gas lifetimes 
should be short. However, at this time, it is difficult to prove 
that the observed C n is not simply the tenuous end product of 
dissociated primordial gas, although the lack of Oi emission 
would be puzzling in this scenario. An alternative scenario 
would be gas production by secondary mechanisms such as 
planetesimal collisions or outgassing from comet-like bodies. 

To calculate the total amount of carbon gas in the disk, 
we assume the disk is similar to the well studied Pictoris 
debris disk. Several gas species have been observed in Pic, 
and their abundance ratios are summarized in Roberge et al. 
(2006). The ratio of neutral and singly ionized carbon in ft Pic is 
Ci/Cii ~ 1. Assuming the same ionization fraction, the total 
column density of carbon in HD 32297 is Nc >5x 10 11 cm -2 . 
The ionization fraction of sodium in Pic is ^0.999 (Roberge 
et al. 2006). This implies a total sodium column density of 
ANa > 2.5 x 10 14 cm' 2 , 500 times larger than the lower limit 
on carbon. 

By assuming solar abundances, we naively expect there to 
be about two orders of magnitude more carbon than sodium. 
This would only be the case if the excitation temperature is 
less than 10 K. We consider three possible explanations. The 
first possibility is that the excitation temperature really is less 
than 10 K; this would be several times lower than the excitation 
temperature measured in Pic (Roberge et al. 2006). Second, the 
disk does not have solar abundances as assumed, but a different 
ratio, meaning less carbon or more sodium. Yet, we expect the 
opposite to the true. Carbon does not feel as strong a radiation 
pressure as sodium because, unlike sodium, carbon does not 
have strong absorption lines in the optical, but in the far-UV 
where the star is much fainter (Roberge et al. 2006). Hence, we 
expect there to be more carbon than sodium relative to solar 
abundances, making the problem worse. The last possibility 
is that ISM sodium absorption lines along the line of sight 
contaminated the HD 32297 sodium measurements, boosting 
the signal. We deem this to be the most likely scenario. 

We can estimate a lower limit on the total gas mass in 
HD 32297 by making a few assumptions. We start by modi- 
fying Equation (1) to get a Cii mass lower limit of Mq u > 
1.7 x 10' 4 M®. By making the same assumptions as above 
(Ci/Cii ~ 1, solar abundances), the total gas mass is ~700 
times the C ii mass, giving a lower limit of M > 0.1 M®. 

Gas in disks can affect the distribution of dust. Gas orbits 
the star at sub-Keplerian speeds due to either to a gas pressure 
gradient or radiation pressure. The dust grains, if large enough 
not to experience strong radiation pressure, orbit at Keplerian 
speeds, and therefore feel a headwind that causes them to spiral 
inward. This mechanism could be a way of transporting dust 
grains from the outer disk to the inner disk. Krivov et al. (2009) 
investigated how 0.3-30 M® of gas would affect HD 32297’s 
radial surface brightness profiles, and found the gas has little 
effect on the disk outside 110 AU. The distribution of the 
inner disk, however, could be affected by small amounts of 
gas. Since we have no limiting spatial information on the gas in 
HD 32297, we cannot determine if the gas significantly affects 
the distribution of the dust in the inner disk. This scenario will 
be further discussed in Section 5.4. 

5.4. Inner Disk 

Our best model for the inner disk of HD 32297 starts at 
~1 AU. The outer edge of the inner disk is unconstrained, but 
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we find a lower limit on the outer edge of ~5 AU. We assumed 
the inner disk grains are astrosilicate grains, and our best fitting 
model has a minimum grain size of 2 jim. The fit to the inner disk 
depends strongly on the fit to the outer disk. Although the results 
of the grain size in the inner disk varied by about two orders 
of magnitude (2.2 /xm-200 // m ) the inner radius only varied by 
a factor of a few (1. 1-3.6 AU). This range places the inner 
disk near the habitable zone of the star. A simple scaling 
of the solar system’s habitable zone (~0.7-1.5 AU) places the 
HD 32297 habitable zone at ~ 1.7-3. 5 AU. 

The presence of dust in the habitable zone does not rule out 
terrestrial planets in this region. Low-mass planets may not have 
had time to clear this region of planetesimals. In fact, some of 
the dust may be trapped in resonance with a planet (Wyatt 2003; 
Stark & Kuchner 2008). 

Since the inner disk has not been resolved, we do not know 
how the dust is distributed. The dust distribution in the inner 
region depends on the location of the parent material and dust 
transportation. We consider two scenarios for the origin of the 
inner dust disk: one where the inner dust disk is fed by dust 
transported from the outer disk through gas drag, and the second 
where there is another planetesimal belt closer to the star. 

In the first scenario, the presence of gas might affect the dust 
distribution, such that it is no longer collision-dominated. The 
entire inner disk may be composed of material that has leaked 
inward from the outer disk due to gas drag. If this were the case, 
we would expect a smooth surface density distribution from the 
outer disk to the inner disk. At first glance, this might seem 
inconsistent with an SED that is well fit by a two-temperature 
blackbody model. This model is most easily interpreted as two 
rings with a gap in between. However, a gap is not needed to 
produce such an SED. The region closest to the star will be 
much warmer and thus brighter than the dust in the intermediate 
region. The signature of dust in the intermediate region would 
be hard to detect in the SED alone. Reidemeister et al. (2011) 
have shown that in the Eps Eri disk, even a bimodal SED 
curve can be reproduced with models that assume transport 
of dust from the outer disk (in that case, caused by stellar 
winds rather than gas), and thus a continuous distribution of 
dust from the outer to the inner region. Since the mass of the 
inner disk that could account for the observed warm emission 
is > 6 x 1 0 9 M®, and assuming the radius of the inner 

disk of ~ 1 AU, such a continuous distribution within ~ 1 00 AU 
would imply roughly 6 x 1CU 5 M® worth of dust. For gas drag 
to work, the gas mass should exceed the dust mass. Therefore, 

0.1 M® of gas, which is a lower limit that we placed from the 
[Cii] observations in Section 5.4, would be sufficient for the 
transport. However, without images of the inner ~50 AU of the 
disk, we cannot confirm gas drag as the origin of the inner dust. 

In our second scenario for the origin of the inner dust, there 
could be another belt of planetesimals closer to the star, similar 
to the asteroid belt. Collisions between planestimals, in this 
case, would produce the dust seen in the SED, just like in the 
outer belt. We made calculations with the model of Lohne et al. 
(2008) and a velocity-dependent critical fragmentation energy 
from Stewart & Leinhardt (2009). They suggest that an asteroid 
belt of a sub-lunar mass at 1.1 AU in a 30 Myr old system 
can easily sustain 1 0 7 A7® to 10 S /V/® worth of dust through a 
steady-state collisional cascade. This is more than the mass of 
the inner disk (M >6x 10~ 9 M®). Without spatial information 
on the inner disk, we cannot tell if the presence of warm dust is 
due to an asteroid belt or from material leaked inward from the 
outer disk. In any case, the mass of the inner disk is significantly 


less than the outer disk; only its proximity to the central star 
makes it outshine the outer disk in the mid-IR. 

6. SUMMARY 

We present new Herschel PACS and SPIRE photometry and 
spectroscopy of the edge-on debris disk around HD 32297 . Our 
main conclusions are the following. 

1 . We detected the disk at 13 wavelengths from 63 to 500 /rm, 
filling in a gap in the SED in this region. The new data 
probe the peak of the thermal emission. 

2. We detected a 3.7er [C n] line at 158 /im, making HD 32297 
only the fourth debris disk with atomic gas detected with 
Herschel. We estimate a lower limit on column density of 
Ncn > 2.5 x 10 11 cm -2 . 

3. The stellar fit to the optical, near-IR and UV data suggest 
the star has a later spectral type than typically quoted, likely 
an Al. 

4. Our SED models require a warm component to fit the large 
mid-IR excess. This material is too warm to be part of the 
ring imaged in scattered light. The geometry of the warm 
component is unconstrained; we were able to fit it with a 
low-density disk at radii greater than 1 AU. 

5. The best fitting model to the outer disk includes grains 
consisting of silicates, carbonaceous material, water ice, 
and a highly porous structure. These grains are similar to 
cometary grains found in the solar system. 

This work is based on observations made with Herschel, a 
European Space Agency Cornerstone Mission with significant 
participation by NASA. Support for this work was provided 
by NASA through an award issued by JPL/Caltech. J.L. and 
J.C.A. thank the PNP-CNES for financial support. We also wish 
to thank the anonymous referee for a very thorough review that 
helped improve this paper. 
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